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ABSTRACT 

New BVRI photometry and optical spectroscopy of the Type Up supernova 2004dj 
in NGC 2403, obtained during the first year since discovery, are presented. The pro- 
genitor cluster, Sandage 96, is also detected on pre-explosion frames. The light curve 
indicates that the explosion occured about 30 days before discovery, and the plateau 
phase lasted about +110 ± 20 days after that. The plateau-phase spectra have been 
modelled with the SYNOW spectral synthesis code using H, Na I, Ti II, Sc II, Fe II 
and Ba II lines. The SN distance is inferred from the Expanding Photosphere Method 
(EPM) and the Standard Candle Method (SCM) applicable for SNe Up. They re- 
sulted in distances that are consistent with each other as well as earlier Cepheid- and 
Tully-Fisher distances. The average distance, D = 3.47 ± 0.29 Mpc is proposed for 
SN 2004dj and NGC 2403. The nickel mass produced by the explosion is estimated 
as ~ 0.02 ± 0.01 Mq- The SED of the progenitor cluster is reanalysed by fitting pop- 
ulation synthesis models to our observed BVRI data supplemented by U and JKH 
magnitudes from the literature. The x^-minimization revealed a possible "young" so- 
lution with cluster age Td = 8 Myr, and an "old" solution with T^i =20-30 Myr. 
The "young" solution would imply a progenitor mass M > 20 M©, which is higher 
than the previously detected progenitor masses for Type II SNe. 
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1 INTRODUCTION 

Type II supernovae (SNe) emerge from shock-driven explo- 
sion of massive (M > 8 Mp^) stars initiated by core collapse 
JWooslev fc WeavedlTosl lNadvozhinll2003l) . In particular, 
plateau Type II SNe (SNe Up) result from core collapse of 
supergiants that have massive hydrogen-rich envelope. The 
light curves of SNe Up are characterized by a constant lu- 
mino sity plateau lasti ng about 80 - 120 days after explo- 
sion JNadvozhinl l2003r . Recently, the progenitors of a few 
SNe Up have been directly identified on pre-explosion im- 
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ages ("Maund fc Smartt"200^: Van Dvk et alJbOOd iLi et alJ 
^005a; Maund. Smartt & Danzigcr 2005; Li et alJ l2005bl) . 
revealing that these are 8 — 15 Mq supergiants, close to 
the low-mass theoretical limit of the core collapse process. 

SN 2004dj in NGC 2403 is an outstanding SN discov- 
ered by Itagaki (Nakano ct al. 2004) on July 31, 2004 (JD 
2453218), because it is one of the brightest and closest su- 
pernova observed ever. Based on its first optical spectrum 
ilPatat et all l2004l) SN 2004dj was quickly classified as a 
"normal" Type Up event. Unfortunately, it was discovered 
during the plateau phase (~ 1 month after explosion), thus, 
the maximum of the light curve co uld not be covered. Du e 
to its proximity (D ~ 3.1 Mpc) (|FVeedniM^^^Jj20M), 
SN 2004di was also detected in radio llStockdale_et_al |20& 
iBeswick et "811120051) . infrared JSueerman fc Van DvlJ Im 
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Figure 1. Left: The field of SN 2004dj observed with the Princeton VersArray 1300B CCD on the Im RCC telescope at Konkoly Obs. 
The field of view shown is about 7x5 arcmins. North is up and East is to the left. The local comparison stars and the SN are labelled. 
The standard magnitudes of the comparison stars are listed in Table 2. Right: The detection of the progenitor cluster Sandage 96 on 
archival frames taken with the 60/90 cm Schmidt telescope and Photometries CCD at Konkoly Obs. on 11th Jan. 2002. The applied 
filter is indicated on the upper left corner on each frame. 



iKotak et al.120 05') and X-ray bands fPool ev fc Lewinl2nn4l) . 
J?-band optical photometry and high-resolution Ha spec- 
tro scopy obtained durinR t he plateau phase was published 
by iKorcakova et alJ l)2n0,'tl . IChugai et alJ (|200,'J l presented 
BVR light curve and full optical spectra in the nebular 
phase. They estimated the amount of ^^Ni as ~ 0.02 Mq, 
and pointed out an asymmetric, bipolar distribution of the 
ejected nickel from modeling the nebular Ha profile. The op- 
tical light curve observed in the BATC (intermediate-band) 
system was also studied by Zhang et al. (2005). They con- 
cluded that SN 2004dj was a typical SN Up that ejected a 
~ 10 Mq envelope and synthesized ^ 0.02 Mq of radioactive 
nickel. 

Particular interest was paid on the association of 
SN 2004dj with the compact cluster Sandage 96 which 
is thought to be th e host of the pro g enitor star 
(lYamaoka et al." "2004"). "Maiz-ApcUaniz e t aP (|20^ and 
[Wang et al. ( 2005) analysed the spectral energy distribution 
(SED) of Sandage 96 (based on archival CCD photometry), 
and concluded that S96 is indeed a compact cluster with age 
of ~ 14 - 20 Myr and total stellar mass of - 10'' - lO'^ Mq 
(the uncertainty is due to the different reddening determi- 
nation). From population synthesis, the initial mass of the 
progenitor star of SN 2004dj was estimated as 12 — 15 Mq, 
which is in the upper part of the mass range of other SNe 
Up progenitors (see above). 

In this paper we present time-resolved optical photome- 
try and spectroscopy of SN 2004dj covering its first year after 
discovery. Our primary goal is to derive a reasonable dis- 
tance to this SN and its host galaxy NGC 2403 by using all 
available information. The new observational data are pre- 
sented in Section 2 and 3. In Section 4 we apply a variant of 
the Expanding Photosphere Method (EPM) and the SNe Up 
Standard Candle Method (SCM) for distance determination, 
and critically compare the results with other distance esti- 
mates of the host galaxy. Based on the improved distance, 
we estimate some of the physical properties of SN 2004dj and 
S96 in Section 5. Finally, Section 6 summarizes our results. 



2 OBSERVATIONS 

2.1 Photometry 

The photometric data were collected with four telescopes, 
whose basic parameters are summarized in Table 1. All ob- 
servations were made through Johnson- Cousins BVRI fil- 
ters. The brightness of the SN was tied to local compari- 
son stars, whose magnitudes were calibrated via Stetson's 
standard field of NGC 7790, when the sky was photometric. 
These local comparison stars are shown in Fig. 1. 

The reduction of the CCD frames was performed in 
IRAF following the standard procedure: bias/dark subtrac- 
tion and fiat field division. The photometry of the reduced 
frames was computed by two methods. First, the instrumen- 
tal magnitudes were derived via aperture photometry. The 
aperture radius was set as 2xFWHM of the stellar images. 
The local background was determined in an annulus with 
inner radius and width of 3xFWHM and 2xFWHM, re- 
spectively. Because of the galaxy contamination on the SN 
and several other local calibrator stars, the background fiux 
around each star has been determined with the built-in cen- 
troid algorithm that is expected to give the best estimate on 
quickly variable background. In order to test this result, the 
background level was also calculated with the mode algo- 
rithm (3 X median - 2 x mean), after eliminating the outliers 
that deviated more than 3(t from the mean value. No signif- 
icant difference was found between the two background lev- 
els, suggesting that the background determination around 
the selected objects is robust. 

The stellar magnitudes coming from the aperture pho- 
tometry were checked with the results of simultaneous PSF- 
photometry. This was done only for the Konkoly Schmidt 
telescope frames, where the larger field of view contained 25 
- 30 bright stars, thus, reliable PSF could be determined. 
This comparison is shown in Fig. 2, where the aperture mi- 
nus PSF magnitudes are plotted against the V magnitude. 
It can be seen that most of the magnitudes of the two pho- 
tometric methods agree well within ±0.05 mag. A few out- 
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Table 1. Basic data of telescopes used for photometric observations. The columns contain the foUowings: code, observatory, diameter of 
telescope, manufacturer of CCD, size of CCD in pixels, readout noise in electrons, pixel scale in arcseconds, field of view on the CCD in 
arcniinutes and typical FWHM of stellar PSFs in arcseconds. Note that the CCD at Szeged Observatory was used with 2x2 binning. 



Code Observatory 



Telescope 



CCD 



Pixels RON (e-) Scale (") FOV {') FWHM (") 



A 
B 
C 
D 



<onkoly 


60/90 cm Schmidt 


Photometries 


1024x1536 


14 


<onkoly 


100 cm RCC 


Princeton VersArray 


1340x1300 


10? 


Szeged 


40 cm Cass 


SBIG ST-9E 


512x512 


13 


FLWO 


120 cm 


Minicam 


4800x4800 


7.6 



1.00 


17x25.5 


3.3 


0.30 


6.7x6.5 


3.0 


0.70 


6x6 


3.8 


0.30 


24x24 


1.2 



0.25 

0.2 

0.15 

0.1 

0.05 



-0.05 

-0.1 

-0.15 
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Figure 2. Difference of the magnitudes obtained from aperture- 
and PSF-photometry against the V magnitude. The symbols de- 
note the filters as follows: filled circles - B, filled squares - V, open 
circles - R, open triangles - I filter. 



Table 2. Standard magnitudes of local comparison stars. Er- 
rors are given in parentheses. See Fig. 1 for object identification. 
The last row lists our measurements of the progenitor cluster 
Sandage 96. 



Star 


V 


B-V 


V -R 


V - I 


A 


16.60(0.05) 


0.71(0.10) 


0.42(0.06) 


0.77(0.05) 


B 


16.80(0.07) 


1.57(0.09) 


0.97(0.05) 


1.84(0.05) 


C 


16.34(0.06) 


0.89(0.07) 


0.57(0.04) 


0.97(0.03) 


D 


14.81(0.04) 


0.76(0.02) 


0.45(0.02) 


0.78(0.02) 


F 


15.61(0.05) 


0.95(0.05) 


0.58(0.05) 


0.97(0.03) 


G 


16.42(0.05) 


0.90(0.06) 


0.53(0.04) 


0.93(0.04) 


H 


12.61(0.03) 


0.98(0.04) 


0.56(0.03) 


1.00(0.02) 


I 


10.06(0.03) 


0.56(0.02) 


0.38(0.02) 


0.64(0.02) 


J 


14.46(0.04) 


0.70(0.04) 


0.43(0.03) 


0.78(0.02) 


K 


14.20(0.04) 


0.74(0.04) 


0.45(0.03) 


0.78(0.02) 


L 


16.66(0.05) 


0.88(0.09) 


0.54(0.06) 


0.92(0.06) 


S96 


17.85(0.05) 


0.40(0.07) 


0.33(0.06) 


0.79(0.05) 



liers correspond to the faintest stars in the B band. It is 
concluded that the results of the aperture photometry with 
the parameters given above are in accord with the PSF- 
magnitudes, therefore, the systematic error caused by the 
photometric method is negligible with respect to the pho- 
ton/detector noise and errors of the standard transforma- 
tion. 

The host galaxy NGC 2403 was also observed in BVRI 
bands with the Schmidt telescope at Konkoly Observatory 
during the course of our nearby galaxy survey programme^ 
(Sarneczkv et al. 2005) on 23 December 2001 and 11 Jan- 
uary 2002. The progenitor cluster Sandage 96 is clearly 
present on these frames (Fig. 1 right panel). Using the same 
set of comparison stars, the magnitudes of S96 were com- 
puted in the same way as described above. The results are 
in the last row of Table 2. Our magnit udes are in very good 
agreem ent with those presented by iMaiz-Apellaniz et alJ 
i2004l) . Because most of the hght of the progenitor cluster 
is still present after the SN explosion, the measured mag- 
nitudes of SN 2004dj have been corrected for the hght of 
the progenitor in the following way: the magnitudes of the 
progenitor as well as those for the SN plus progenitor have 
been converted into fluxes, then the progenitor fluxes have 
been subtracted from the SN-|-progenitor fluxes. The result- 
ing fluxes, expected to be due to only the SN, were trans- 
formed back into magnitudes. 

In order to further test the reliability of the back- 
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Figure 3. Results of the image subtraction method illustrated 
with the R frames obtained on May 5, 2005. Left panel: the broad- 
ened object frame after convolving with the kernel produced by 
psf match. Middle panel: the template frame scaled to the light 
level of the object frame. Right panel: the difference between the 
two images. The galaxy and the unsaturated stars have been re- 
moved. 



ground subtraction and the correction for the progenitor, 
the Konkoly Schmidt frames were also analysed using the 
image subtraction method. This was also done with built-in 
IRAF tasks. Three epochs have been selected: Aug. 8, 2004 
(-f39 days after explosion), Nov.02, 2004 (-1-125 days) and 
May 13, 2005 (-1-317 days), when the SN was in the middle 
of the plateau phase, during the transition phase and in the 
nebular phase, respectively. The template frames containing 
the progenitor have been subtracted from the BVRI frames 
obtained on these epochs with the Schmidt telescope in the 
following way. First, the images have been registered with 
the geomap, geotran and imalign tasks. Second, the stellar 
PSFs on the object and template frames have been matched 
using the psf match task. Because the PSFs of the template 
frames were broader than those of the object frames, the lat- 
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Figure 4. Magnitude corrections due to the progenitor cluster 
computed with i) direct aperture photometry of the SN and tem- 
plate frames (filled circles) and ii) image subtraction method 
(open crcles) as a function of wavelength. The phase of the obser- 
vation is labelled on the right. The magnitude corrections from 
the two methods agree within 0.02 mag, even at the last epoch, 
when the SN was the faintest on our frames. 



ter images have been broadened with the kernel produced 
by psf match. Third, the intensities of the comparison stars 
have been scaled to the same level on both frames with the 
task linmatch. Fourth, the scaled template frames have been 
removed from the broadened object frames to produce the 
subtracted images. The result of this procedure is illustrated 
in Fig. 3. It is seen that the galaxy and the unsaturated stars 
have been nicely removed by the image subtraction, leaving 
a well-defined object (the SN) on a flat background. Then, 
the aperture photometry have been computed again, now 
measuring the comparison stars on the broadened object 
frames, and the SN on the subtracted frames. The SN mag- 
nitudes obtained this way were compared with the results 
of the first method, i.e. that were computed directly on the 
object frames and corrected for the progenitor magnitudes 
(see above). Fig. 4 shows the difference between the uncor- 
rected and progenitor-corrected SN magnitudes in the case 
of direct aperture photometry (filled circles), and the dif- 
ference between the SN magnitudes measured before and 
after the template subtraction (open circles), as a function 
of wavelength, for the three epochs considered. It is visible 
that the corrections coming from the direct photometry do 
agree within 0.02 mag with those from the image subtraction 
method, even at the last epoch when the SN was the faintest 
on our frames. This also means that the galaxy contam- 
ination could also be effectively removed with the selected 
aperture/annulus combination in the first method. Note that 
this good agreement is partly due to the fact that the object 
and template frames have been made with the same instru- 
ment. The image subtraction method is expected to perform 
less well for those SN frames that have been made with 
other telescopes. Moreover, the digital image subtraction al- 
ways increases the noise on the difference image, resulting in 
slightly higher uncertainties. Therefore, it is concluded that 
the light of the progenitor cluster can be reliably removed 
from the SN magnitudes obtained with direct aperture pho- 
tometry of the object frames (containing the galaxy and the 
SN) using the progenitor magnitudes listed in Tabic 2. The 
image subtraction method (that removes the host galaxy 
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Figure 5. Top panel: Light curve of SN 2004dj in BVRI bands. 
The curves have been shifted vertically for better visibility with 
the amount indicated in the legend on the right. Bottom panel: 
The light curve during the early tail phase. The line indicates the 
expected slope of the Co-decay. 



and the progenitor pixel-by-pixel) gave essentially the same 
result for the selected frames. In the foUowings we use the 
magnitudes of SN 2004dj computed with the first method 
and corrected for the progenitor magnitudes in Table 2. 

Finally, all photometric data were transformed into the 
standard Johnson- Cousins system using colour terms de- 
termined for each telescope/CCD on photometric nights. 
The final calibrated and progenitor-corrected magnitudes of 
SN 2004dj are listed in Table 3. The errors (given in paren- 
theses) reflect the observational noise and the uncertainties 
of the magnitudes of local comparison stars (see Table 2) . 

The light curves of SN 2004dj are plotted in Fig. 5. This 
is a typical Type Up light curve. During the plateau phase 
the light level is almost constant in V, R, I bands. There is 
a slow decrease in B, which would be even stronger in the 
U band. It is generally accepted that this plateau phase is 
due to a hydrogen recombination front that moves inward 
the expanding H-rich ejecta. Then, following a rapid transi- 
tion, the SN enters the nebular phase (often referred as the 
tail phase), where the light variation reflects the radioac- 
tive decay of ^^Co into ®^Fe. The bottom panel of Fig. 5 
shows only the radioactive tail with the expected slope of 
the Co-decay. It is visible that in the early tail phase the 
slope of the light curve is less than that of the Co-decay in 
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Table 3. Observed magnitudes of SN 2004dj. The columns contain the followings: date and JD-2450000 of the observations, days relative 
to the estimated moment of explosion (JD 2453187, see Sect. 4), BVRI magnitudes corrected for the magnitudes of S96 (errors are in 
parentheses) and codes of the applied telescopes. See Table 1 for the telescope codes. 



Date 


JD-2450000 


* — texpl 


B 


V 


R 


/ 


Tel. 


2004-08-05 


3223.4 


36 


12.70 


0.05) 


12.05 


0.04) 


11.61 


0.03) 


11.50 


0.07) 


C 


2004-08-08 


3226.6 


39 


12.73 


0.06) 


12.00 


0.03) 


11.57 


0.03) 


11.40 


0.03) 


A 


2004-08-10 


3228.4 


41 


- 




12.02 


0.03) 


11.57 


0.03) 


11.45 


0.05) 


C 


2004-08-10 


3228.6 


41 


12.81 


0.05) 


12.00 


0.03) 


11.58 


0.03) 


11.39 


0.03) 


A 


2004-08-11 


3229.6 


42 


12.81 


0.05) 


12.01 


0.03) 


11.58 


0.03) 


11.39 


0.03) 


A 


2004-08-12 


3230.4 


43 


12.69 


0.05) 


11.99 


0.03) 


11.54 


0.05) 


11.39 


0.03) 


C 


2004-08-16 


3234.4 


47 


13.06 


0.07) 


12.04 


0.05) 


11.61 


0.03) 


11.42 


0.03) 


C 


2004-08-18 


3236.6 


49 


12.98 


0.06) 


12.04 


0.03) 


11.61 


0.03) 


11.40 


0.03) 


B 


2004-08-28 


3246.4 


59 


12.89 


0.12) 


12.10 


0.03) 


11.64 


0.03) 


11.43 


0.03) 


C 


2004-08-29 


3247.3 


60 


13.07 


0.06) 


12.11 


0.03) 


11.68 


0.03) 


11.45 


0.03) 


C 


2004-09-17 


3266.7 


79 


13.38 


0.05) 


12.33 


0.03) 


11.88 


0.03) 


11.65 


0.03) 


A 


2004-11-02 


3312.4 


125 


16.16 


0.05) 


14.93 


0.04) 


14.36 


0.03) 


13.95 


0.03) 


A 


2004-11-11 


3321.4 


134 


16.28 


0.05) 


15.01 


0.04) 


14.39 


0.03) 


13.99 


0.03) 


A 


2004-11-15 


3325.7 


138 


16.31 


0.05) 


15.04 


0.04) 


14.43 


0.03) 


14.03 


0.03) 


A 


2004-11-24 


3334.6 


147 


16.41 


0.17) 


15.15 


0.08) 


14.49 


0.08) 


14.12 


0.06) 


B 


2004-12-04 


3343.5 


156 


16.47 


0.05) 


15.21 


0.04) 


14.48 


0.04) 


14.11 


0.03) 


B 


2004-12-07 


3346.7 


159 


16.47 


0.05) 


15.24 


0.04) 


14.46 


0.04) 


14.10 


0.03) 


B 


2004-12-29 


3369.6 


182 


- 




15.46 


0.06) 


14.58 


0.04) 


14.18 


0.03) 


A 


2004-12-30 


3370.6 


183 


16.64 


0.05) 


15.47 


0.04) 


14.58 


0.05) 


14.17 


0.03) 


A 


2005-01-05 


3376.4 


189 


16.61 


0.05) 


15.50 


0.04) 


14.60 


0.05) 


14.17 


0.03) 


A 


2005-01-21 


3391.5 


204 


16.77 


0.05) 


15.66 


0.04) 


14.68 


0.04) 


14.28 


0.03) 


A 


2005-02-10 


3412.4 


225 


- 




15.86 


0.06) 


14.83 


0.04) 


14.43 


0.05) 


A 


2005-02-26 


3428.4 


241 


16.90 


0.06) 


15.93 


0.04) 


14.88 


0.05) 


14.45 


0.05) 


A 


2005-02-28 


3430.4 


243 


16.99 


0.06) 


15.90 


0.04) 


14.89 


0.04) 


14.51 


0.03) 


A 


2005-03-31 


3461.3 


274 


17.10 


0.07) 


16.19 


0.05) 


15.12 


0.04) 


14.74 


0.04) 


A 


2005-04-02 


3463.3 


276 


17.11 


0.07) 


16.16 


0.05) 


15.12 


0.04) 


14.75 


0.04) 


A 


2005-04-04 


3465.3 


278 


17.21 


0.07) 


16.21 


0.05) 


15.13 


0.04) 


14.76 


0.04) 


A 


2005-04-05 


3466.3 


279 


17.13 


0.08) 


16.24 


0.05) 


15.17 


0.04) 


14.82 


0.04) 


A 


2005-05-12 


3502.1 


315 


17.44 


0.10) 


16.66 


0.05) 


15.55 


0.04) 


15.02 


0.04) 


D 


2005-05-13 


3504.3 


317 


17.36 


0.10) 


16.46 


0.05) 


15.44 


0.04) 


15.10 


0.04) 


A 



R and /. The decline rate in these bands is 0.6 niag/(100 
days) compared with the 0.98 mag/(100 days) rate of the Co- 
decay. On the other hand, the B and V hght curve declines 
with 0.8 mag/(100 days) which is c lose to the expected rate. 
Similar phenomenon was found bv lZhang et alJ l|2005|) from 
their intermediate-band hght curves: in the AA6600 — 8500 
A wavelength regime the decline rates are much less than 
the Co-rate. Since the R and / bands are dominated by the 
Ha and Ca-triplet emission lines, which can be particularly 
strong at the beginning of the tail phase, it is probable that 
the reduced decline rate in these bands reflects the contribu- 
tion of these lines to the thermal radiation powered by the 
trapping of gamma-rays and positrons from the Co-decay. 



IChueai et alJ (|2005|) also presented a tail light curve 
of SN 2004dj obtained through Johnson- Cousins B, V and 
R filters. Comparing their published photometry with ours 
presented here, the hght curve slopes in different bands are 
consistent. However, there is a general, systematic offset of 
0.3 - 0.4 mag in all bands in the sense that the Chugai et al. 
magnitudes are systematically brighter. Both sets of hght 
curves have been transformed to the standard system, al- 
though there is no st ar in common in th e two groups of 
local calibrators. Since IChueai et al.l ll2005h do not mention 
whether they had corrected their data for the progenitor 
light, it is probable that this difference is partly due to the 
presence of the progenitor hght in their magnitudes. But it 



is interesting that if we do not correct our data for the pro- 
genitor, the offset is still present, although with a reduced 
amount (Chugai et al.'s data are still brighter by ~ 0.2 mag). 
At present, it is not clear what is the cause of this systematic 
offset, but note that the standard transformation of SN mag- 
nitudes obtained with different telescope/CCD/filter com- 
binations are always uncertain at the level of ~ 0.1 mag 
JSuntzefJl2OO0l) due to the nonstellar SED of the SN, espe- 
cially during the nebular phase. This effect may easily be 
responsible for most of the remaining offset. 



2.1.1 Comparison with other SNe Type lip 

ft is interesting to compare the hght curve of SN 2004dj 
with those of other Type Up SNe, although it is known 
that SNe U p light curves are quite hctcrogcnious. Recently 
iHamuv fc P into (2002) and ffamuv (2005) showed that 
there is a correlation between the absolute magnitudes in 
the middle of the plateau and the photospheric radial veloc- 
ity at the same epoch, which partly explains the observed 
difference s in the light curves of different SNe ffp. 

Chu gai et alj (2005) compared SN 2004dj with 
SN 1999gi, another well-observed Type Up SN. Based on the 
V Hght curve (supplemented by amateur observations) and 
an estimated explosion date JD 2453170, th ey concluded 
that the two light curves are very similar. IZhang et alJ 
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Table 4. Light curve parameters of Type lip SNc used for com- 
parison. 
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Figure 6. Comparison of the hght curves of SN 2004dj (filled 
circles), SN 1999gi (open circles), SN 1999em (open squares) and 
SN 2003gd (asterisks). The magnitudes have been corrected for 
reddening and distance listed in Table 4. The curves have been 
shifted vertically for better visibility (the amount of the applied 
magnitude shift is indicated on the left-hand side of each curve 



IJ200)tI) did the same but using the data of SN 1999ein for 
comparison. Based on the assumption that the two hght 
curves should be similar, they estimated the moment of ex- 
plosion as JD 2453 167 ± 21, which was nearly the same as 
IChugai et al.l (|200.'j) determined. 

We compared our new light curve (Fig. 5 and Ta- 
ble 3) with the pu blished magnitude s of SN 1999em 
l|Leonard et alJ I2 n02al) . SN 1999gi (Leonard et al.l l2002bl ) 
and SN 2003gd JHendrv et al.l |2005l) To do this, aU data 
were corrected for reddening and distance modulus, col- 
lected from references listed in Table 4. F or SN 1999em we 
have us ed the distance derived recently bv lPessart fc Hillien 
I|2n05a) . which is in good agreement with the Cepheid- 
based distance of its host galaxy. For SN 2004dj we have 
applied our preferred values of E[B — V) = 0.07 mag, 
Texpi = 2453187 and D = 3.47 Mpc (see Section 4). The ab- 
solute magnitudes for each SN are plotted together in Fig. 6. 

From Fig. 6 it is visible that the absolute light level 
of SN 2004dj is similar to that of SN 1999gi during the 
plateau phase, but different from those of SNe 1999em and 
2003gd. This difference may be due to different nickel masses 
for these SNe. However, as it will be discussed in Section 
2.2, the spectra and the expansion velocities are quite sim- 
ilar for SN 2004dj and SN 1999em in the middle of the 
platea u phase, tha t would suggest similar plateau magni- 
tudes JHamuv fc Pinto 2002) . It seems that the plateau 
phase of SN 2004dj ended slightly early relative to 99em 
and 99^1. Com paring our data with the 7? light curve of 
iKorcakova et al. (2005), we estimate the starting date of 
the transition into the nebular phase as i = + 80 ± 20 
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days (the uncertainty is due to both the poor sampling of 
the light curves around this epoch and t he uncertainty of 
the explosion date). The inflection time IIElmhamdi et all 
■2003a1 for SN 2004dj is estimated as ti = -1-110 ± 20 days, 
which is shghtly earher than that for SN 1999em (+120 ± 4 
days) and SN 1999gi (-1-125 ± 3 days), but this is not a sig- 
nificant difference because of the relatively high uncertainty 
of the explosion date of SN 2004dj. 

In the nebular phase, SN 2004dj has the faintest light 
level among these four SNe, except in B where SN 2003gd 
is slightly fainter. The faintness of SN 2004dj is increasing 
toward longer waveleng ths. In this phase SN 2004di is more 
simil ar to SN 2003gd (IVan Dvk et alJ 120031: [Hendry et alJ 
I2OO5I) than the other two SNe. lElmhamdi et al.l tOOSsf 
found that the tail light level of SNe Up correlate well with 
the plateau magnitude in V . From Fig. 6 it is suspected that 
this correlation is actually weak in V, and probably not valid 
at all in R and /. 

Since SNe Up tail light curves are powered by ra- 
diactive decay of ^^Co - ^^Fe, different brightness prob- 
ably reflects the difference in the amount o f the syn- 
thesi z ed iron-peak elements, particularly ^^Ni JNadvozhinl 
|2003|: [Elmh amdi et alJl2003a^ . From hydrodynamic models 
iNadvozhinI (l2003h pointed out that the middle-plateau abso- 
lute magnitude My is proportional to the plateau duration 
Atp, the expansion velocity in the mi ddle of the plateau 
and the explosion energy (see Eq.l of INadvozhinI (12003.1 ). 
Thus, from the shape of the light curve alone, we suspect 
that SN 2004dj had similar explosion energy and nickel mass 
than those of SN 1999gi, but less than those of SN 1999em. 
We will examine this question more quantitatively in Sec- 
tion 5. 



2.1.2 Reddening 

Reddening is the key parameter when photometry is used 
to derive physical parameters. From the literature, the red- 
dening of SN 2004dj is controversial, fatal ct al. (200J) re- 
ported the identiflcation of interstellar Na D (from an un- 
resolved absorption trough) in the very first spectrum of 
SN 2004dj, from which they deduced E(B - V)totai = 0.18 
mag. 1 day later iGuenther fc Klosd (|200^ stated that with 
high-resolution echelle spectroscopy they resolved the Na D 
lines in the host galaxy NGC 2403. They revealed the con- 
tribution of the host galaxy to the reddening of SN 2004dj 
as E{B ~V)host ~ 0.026 mag. Sin ce the galactic comp onent 
is only E(B - V)gai = 0.04 mag iSchlegel et al.lll998l) . this 
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Figure 7. The colour curves of SNe 2004dj, 1999em, 1999gi and 
2003gd corrected for galactic reddening (see text). The symbols 
are the same as in Fig. 6. 



indicates si gnificantly le s s tota l reddening {E{B — V)totai = 
0.066) than JPatat et alJ l|2004l estimated. 

Reddening has also been determined indirectly from 
population synthesis models of the progenitor clus- 
ter S96. Fitting cluster SEDS to archival optical 
and near-IR (UBVRIJH K) photometric data of S96, 
iMaiz-Aoellaniz et alJ l)2004) defwed E(B ~V)totai ~ 0.13± 
0.03 mag. This value is between the two direct spectral 
measurement mentioned abov e, being def i nitely in better 
agreement with the result of IPatat et alJ l|2004h . On the 
other hand. lWang et alJ l|2005|) have applied essentially the 
same technique using their own O-IR photometric data 
(through intermediate-band filters, providing better wave- 
length coverage in the optical) on S96. They have obtained 
E{B — V)totai = 0.35±0.05 mag, arguing that this value pro - 
vides a better fit of the calculated SEDs. lZhane et al.l (|200.4l 
also derived E{B — V)totai = 0.33 ± 0.11 mag from the re- 
an alysis of the EW of the interstellar Na D trough reported 
bv lPatat et al.l i2004h . It is clear that the present estimates 
of the extragalactic reddeni ng toward SN 2004di ran ge from 
an almost negligible value (l&^gnthg^fc Klggal2QQ ^ to sig- 
nifica nt in- host extinction (|Wang et al.M200!TL IZhang et alJ 
|2003). 



11 

12 

„ 13 

CD 
(0 

% 14 

" 15 

16 

17 



SN 2004dj 



50 100 150 200 250 300 
Time since explosion (days) 



350 



Figure 8. UVOIR bolometric magnitudes of SN 2004dj, cor- 
rected for reddening {E(B — V) = 0.07 mag). The bolometric 
fluxes (Table 5) have been transformed into magnitudes via Eq.l. 



Comparing our observed colour curves of SN 2004dj 
with those of SNe 1999em, 1999gi and 2003gd some con- 
straints ca n be determined on th e upper or lower limits of 
reddening llLeonard et alj|2002bl) . Fig. 7 shows the B -V, 
V — R and V — I curves of these four SNe Up corrected for 
the galactic reddening (the following values were adopted 
from NED: E{B - V)gai = 0.04, 0.04, 0.017 and 0.069 mag 
for SNe 2004dj, 1999em, 1999gi and 2003gd, respectively). 
SN 2004dj has the bluest colours among these SNe, espe- 
cially in the transition and the nebular phase. Since SNe 
1999em, 1999gi and 2003gd suffered only moderate redden- 
ing in their hosts (see Table 4) , the bluer colour of SN 2004dj 
means that E{B — V)totai > 0.1 mag is less probable. Thus, 
the observed colour curves favour lower reddening, such as 
that determined by Guonthcr & Klosc (2004). This may also 
be preferred, because this is the only one that is based on 
direct measurement of the resolv ed interstellar Na D lines 
within the host galaxy. Note that IChugai et al.l (I2OO3) also 
selected E{B - V)totai = 0.062 mag for SN 2004dj based on 
their observed colour indices, but they attributed it entirely 
to galactic extinction. 

In the foUowings we adopt E{B — V)totai = 0.07 ± 0.1 
mag as the most probable total (galactic plus host) redden- 
ing of SN 2004dj. Its uncertainty is estimated from the scat- 
tering of the spectroscopic measurements mentioned above. 



2.1.3 The bolometric light curve 

The UVOIR bolometric light curve has been calculated 
from the observed BVRI magnitudes plus extrapolating 
an assumed (over-)simplified SED to the wavelength bands 
not covered by the observations. First, the observed magni- 
tudes were corrected for extinction using E{B — V) = 0.07 
mag (see pr evious section) and the standard galactic red- 
dening law JSchlegel et aljfigg^ . ii"— corrections were ne- 
glected, because the redshift of the host galaxy is very 
small, z = 0.00044 (NED). The dereddened magnitudes were 
then transformed to fluxes using the calibration given by 
iHamuv et alJ i200l|) . The fluxes have been integrated nu- 
merically using the effective wavelengths and FWHM values 
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Table 5. UVOIR bolometric fluxes and magnitudes of SN2004dj. 



(10- 



JD-2450000 


t - texpl 




(days) 


3223.4 


36 


3226.6 


39 


3228.6 


41 


3229.6 


42 


3230.4 


43 


3234.4 


47 


3236.6 


49 


3246.4 


59 


3247.3 


60 


3266.7 


79 


3312.4 


125 


3321.4 


134 


3325.7 


138 


3334.6 


147 


3343.5 


156 


3346.7 


159 


3370.6 


183 


3376.4 


189 


3391.5 


204 


3428.4 


241 


3430.4 


243 


3461.3 


274 


3463.3 


276 


3465.3 


278 


3466.3 


279 


3502.1 


315 


3504.3 


317 



fbol 


mtoi 


10 erg/s/cm2) 


(mag) 


4.14 


11.94 


4.35 


11.89 


4.30 


11.90 


4.30 


11.90 


4.44 


11.87 


4.02 


11.97 


4.12 


11.95 


4.06 


11.97 


3.88 


12.02 


3.16 


12.24 


0.33 


14.68 


0.32 


14.74 


0.31 


14.77 


0.28 


14.86 


0.28 


14.86 


0.28 


14.86 


0.25 


14.97 


0.25 


14.97 


0.23 


15.08 


0.19 


15.26 


0.19 


15.31 


0.15 


15.53 


0.15 


15.53 


0.15 


15.56 


0.14 


15.59 


0.11 


15.86 


0.11 


15.86 



Table 6. Journal of spectroscopic observations. The columns con- 
tain the following data: date and JD (— 2450000) of the obser- 
vations, days relative to the estimated moment of explosion (JD 
2453187, see Sect.4.1), central wavelength and airmass. 



of the BVRI filters by a simple midpoint rule. Tlie fiuxes 
in the missing UV and fR bands were extraplolated linearly 
from the B and I fiuxes assuming zero flux at 3400 A and 
23000 A. 

This approximation has been tested by integrating 
the known BVRI magnitudes of Vega and the Sun 
JHamuv et al.l l2001^ and comparing the resulting quasi- 
bolometric fluxes with the observed ones. The relative error 
of the bolometric flux calculated in this way turned out to be 
9.6 % and 0.6 % (0.11 and 0.001 mag) for the Sun and Vega, 
respectively. Because during the plateau phase the SED of 
the SN is similar to stellar photospheres, and the effective 
temperature in most cases is between the temperature of the 
Sun ad Vega (5800 < Tsn < 10000 K), we believe that 
the uncertainty of our bolometric fluxes due to the approx- 
imation applied does not exceed significantly the level of 10 
percent at the early phases. 

During the nebular phase, the SED of the SN becomes 
nonstellar, therefore the comparison with stellar fiux dis- 
tributions may be misleading. Fortunately, SN 2004dj was 
detected with the IRAC and MIPS instruments onboard the 
Spitzer Space Teles cope on four epoch s at the beginning of 
the nebular phase. iKotak et al.l (|2003) presented the pho- 
tometry made in 5 channels in the mid-IR regime (from 3.6/i 
to 24/i). These infrared fluxes were used to test the quality 
of our quasi-bolometric fluxes in the nebular phase. Assum- 
ing negligible reddening in the mid-IR regime, the SED of 
SN 2004dj w as constructed usin g the IR fluxes obtained on 
Nov.l, 2004 llKotak et alJl2005^ and our optical fluxes ob- 
served on the same epoch (Table 3). The missing near-IR 



Date 


JD-2450000 


Days 


Xc 


Airmass 


2004-08-17 


3234.9 


47 


5800 


1.57 


2004-08-20 


3237.6 


50 


5800 


2.92 


2004-08-22 


3239.9 


52 


5800 


1.47 


2004-08-25 


3242.9 


55 


5800 


1.49 


2004-09-03 


3251.8 


64 


5800 


1.49 


2004-09-04 


3252.8 


65 


5800 


1.53 


2004-09-06 


3254.9 


67 


5800 


1.04 


2004-09-22 


3270.9 


83 


5800 


1.22 


2004-10-02 


3281.7 


94 


5800 


1.24 


2004-10-03 


3282.9 


95 


5800 


2.02 


2004-10-05 


3284.9 


97 


5800 


2.25 


2004-10-06 


3285.5 


98 


5800 


1.13 


2004-10-07 


3286.9 


99 


5800 


1.70 


2004-11-14 


3324.7 


137 


5800 


1.32 



fluxes were estimated by interpolation between the optical 
and mid-IR fluxes, assuming power-law wavelength depen- 
dence. This SED wcis then integrated along wavelength us- 
ing the same approximation in the UV boundary as before, 
but extending the IR boundary up to 8 microns. The result 
was 0.317 • 10"^° erg/s/cm^ comparing with 0.33 • 10"" 
erg/s/cm^ estimated from the optical data. The relative dif- 
ference in only 4 percent, which is nearly the same as it 
was in the photospheric phase. Thus, it is concluded that 
the approximative, quasi-bolometric fluxes, estimated from 
BVRI fluxes have about 10 percent uncertainty, both in the 
plateau and the early nebular phase. 

Finally, the bolometric fluxes have been transformed 
into magnitudes using 



mboi 



2.5 logj,, [fboi_ 



11.512 



(1) 



where the zero-point was determined from the solar lumi- 
nosity and absolute bolometric magnitude {Mi,oiq ~ +4.72 
mag and Lq — 3.847 • 10^^ erg s'^ was used). The bolomet- 
ric magnitudes are plotted against time in Fig. 8. The slope 
of the light curve is clearly increasing in the radioactive tail 
phase (as it was also presented in Fig. 5). At t « -1-120 days 
after explosion the slope is about —0.45 mag/100 days, while 
around t ~ -1-300 days it increases up to —0.8 mag/100 days 
which is closer to the expected theoretical value of the Co- 
decay (-0.98 mag/100 days). We deduce that SN 2004dj 
entered fully the radioactive nebular phase around -1-350 - 
-1-400 days after explosion. We analyse the bolometric light 
curve further in Section 4 and 5. 



2.2 Spectroscopy 

The spectroscopic observations were conducted at David 
Dunlap Observatory (DDO) with the 74" telescope and the 
Cassegrain spectrograph. The 100 lines/mm grating was ap- 
plied giving a 2-pixel-resolution of about 800 at 6000 A. The 
journal of spectroscopic observations is presented in Table 6. 
The spectra were reduced and calibrated in the 
usual way using IRAF. After bias, flatfield and cosmic- 
ray corrections, the spectrum was extracted with the 
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Figure 9. Spectral evolution of SN 2004dj during the plateau 
and the transition phase. Labels indicate the elapsed time since 
the estimated date of explosion in days. The spectra have been 
shifted vertically by 5 ■ 10"^'* erg/s/cm^/A for better visibility. 



noao.twodspec.apall task. The wavelength calibration was 
done using two exposures of a FeNe spectral lamp made 
before and after each target obs ervation. For fl u x cal i- 
bration, HP 428 1 8 (sp type AOV, lAdelman et"!!] ijlQSfltl : 
iKharitonov et alJ il988l) ) was observed as a flux standard. 

The DDO Cassegrain spectrgraph has a known limita- 
tion in obtaining low resolution spectra, because the slit can- 
not be rotated, therefore the observations cannot be made 
at the parallactic angle. Due to this technical limitation, 
the continuum slope of our spectra is certainly affected by 
the differential atmospheric refraction, especially in the blue. 
The most distorted spectrum is the one that was taken on 
20th August, 2004, when the airmass was nearly 3 (see Ta- 
ble 5). On the other hand, those spectra that were observed 
between airmasses of 1 and 1.5 are probably less distorted, 
thus, they may be used for quantitative analysis. Neverthe- 
less, the blue end of all spectra was set as 4300 A in order 
to eliminate the most unreliable blue part. 

In the following subsections we present and discuss the 
plateau and nebular spectra of SN 2004dj separately. 



2.2.1 Plateau phase 

Fig. 9 shows five representative spectra from the 13 ones ob- 
tained during the plateau phase. The plateau-phase spectra 
are typical for Type Up SNe: dominant hydrogen lines with 
P Cygni profiles and moderately strong metallic line blends 
in the blue, superimposed on a continuum that is increasing 
toward the blue. Since the spectral evolution is slow in this 
phase, the other spectra are very similar to those plotted in 
Fig. 9. 




4000 4500 5000 5500 6000 6500 7000 7500 8000 
Rest wavelengtti (A) 

Figure 10. Comparison of the spectra of SN 2004dj (thick line) 
and SN 1999em (thin line). The time since explosion is labelled. 
The spectra were dereddened, corrected for redshift of the host 
galaxy and shifted vertically for better visibility. 



The evolution of the spectral features is also similar to 
other SNe Type Up. When the ejecta cools, the continuum 
gets fainter and redder. The emission component of Ha be- 
comes stronger, while the strength of HP and the other ex- 
cited Balmer lines decreases. Among the metallic lines, Na 
D develops a strong P Cygni profile at the end of the plateau 
phase. Unfortunately, no spectrum earlier than -1-47 day was 
available for us, so we could not check the presence of He in 
the early phases. 

Our first spectrum of SN 2004dj was compared with the 
spectra of SN 1999em jLeonard et al. 2002a) made at similar 
phases. These comparison spectra were downloaded from the 
SUSPECT ^ database. Fig. 10 shows that the spectra of the 
two SNe are very similar. This is slightly in contrast with 
the difference between the light curves of the two SNe found 
in Sect. 2. 1.1. 

From the photospheric spectra the radial velocities of 
the expanding ejecta were estimated fr om the Doppler- 
shift of the A5169 A Fe II line. Recently IPessart fc HiUieil 
(|2005a|) pointed out that this line is a good indicator of 
the photospheric velo city. This line is also favoured by 
iHamuv fc Pintd l|2003) for estimating the expansion velocity 
in the middle of the plateau phase. The velocities are plotted 
in Fig. 11 as filled circles. Because our spectra do not cover 
the plateau phase well, the first radial velocit y point was 
computed from the line positions reported by IPatat et alJ 
12004). Fig. 11 also contains the radial velocities estimated 
from the minimum of Ha. As it is expected, the optically 
thick Ha forms at higher atmospheric levels, probing higher 
ejecta velocities than the photospheric Fe II lines. 

The measured radial velocities must be corrected for the 
motion of the host galaxy and the SN itself within the host 
(the barycentric motion of the Earth is usually negligible in 
the case of SNe). This is often done using the redshift of 
the host galaxy adopted from major databases. In the case 
of SN 2004dj, the host galaxy, NGC 2403, has a recession 
velocity of -1-131 km/s (NED). Moreover, the SN is located 



http://bruford.nhn.ou.edu/~suspect/indexl.html 



10 J. Vinko et al. 



14000 
12000 
10000 

5^ 8000 

E 
:^ 

>- 6000 
4000 
2000 



■ 












■ 
■ 












□l] 












- -» 


Ha 








- 


m 













n 

Fell A5169 




•I 







20 40 60 80 

Time since explosion (days) 



100 



120 



Figure 11. Radial velocities of SN 2004dj. The photosplieric ve- 
locities (filled circles) have been computed from the Doppler-shift 
of the A5169 Fe II line (see text). For c omparison, pho t ospheri c 
velocities of SN 1999em (open squares, [L eonard ct al.' 120023^ ') 
and SN 1999gi (filled squares, iLeonard et al. (2002b) ) are also 
plotted. The Ha velocities of SN 2004dj are shown as open cir- 
cles. 



Table 7. Photospheric radial velocities of SN 2004dj corrected for 
redshift (-fl31 km/s) and host galaxy rotation (+90 km/s). The 
third column lists the average radial velocities calculated from 
the A5169 Fe II line. The fourth column contains the velocities 
based on the minimum of Ha. Errors are given in parentheses. 

JD-2450000 Days Vr (km/s) d/jq (km/s) 



3220.671 


33 


4150 


(300) 


6790 (300) 


3234.89 


47 


3423 


(250) 


5693 (300) 


3237.59 


50 


3236 


(250) 


5382 (300) 


3239.89 


52 


3184 


(250) 


5444 (300) 


3242.89 


55 


3183 


(250) 


5353 (300) 


3251.86 


64 


2888 


(250) 


5001 (300) 


3252.85 


65 


2917 


(250) 


4700 (300) 


3254.86 


67 


2813 


(250) 


4577 (300) 


3270.90 


83 


2592 


(250) 


4589 (300) 


3281.67 


94 


2342 


(250) 


4373 (300) 


3282.87 


95 


2214 


(250) 


4312 (300) 


3284.87 


97 


2239 


(250) 


4428 (300) 


3285.49 


98 


2136 


(250) 


4443 (300) 


3286.86 


99 


2230 


(250) 


4256 (300) 



from lPatat et alj i200' 



close the outer end of a spiral arm, thus, its radial velocity 
is expected to differ from the systemic velocity of the whole 
galaxy. Indeed. iFraternali et alJ (|200J) recently mapped the 
kinematics of NGC 2403 via H I emission. It is visible in their 
Fig. 1 that the area around SN 2004dj is on the receding side 
of the galaxy. The excess velocity was estimated to be -1-90 
km/s, thus, the total radial velocity of the barycenter of 
SN 2004dj was set as -1-221 km/s. The final radial velocities 
of the SN ejecta were calculated with respect to this value 
(see Table 7). 
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Figure 12. A nebular spectrum of SN 2004dj, +137 days after 
the explosion. 



2.2.2 Nebular phase 

Fig. 12 shows a single spectrum obtained at the begin- 
ning of the nebular phase {t = +137 days after explosion). 
The spectrum is dominated by strong Ha emission and P 
Cygni profiles of Na I, O I and Fe II. The for bidden [Ca II] 
line is also apparent around 7300 A. Recently IChueai et alJ 
(J200a) analysed a more extensive set of nebular spectra of 
SN 2004dj. They have pointed out significant asymmetry of 
the Ha profile, from which they concluded that the ejec- 
tion of ^®Ni (which is the primary source of Ha emission in 
the nebular phase via radioactive decay) was bipolar inside 
a spherically symmetric envelope. We do not consider this 
topic in detail, except noting that our nebular spectrum of 
SN 2004dj (Fig. 12) do confirm the strong bl ueward asym- 
metry of the Ha emission profile found by IChueai et alJ 
l|2005h . So far this phenomenon was reported only i n two 
other SNe, nam ely SN 1987A Jphillips fc Williams! Il99 J) 
and SN 1999em jEknhamdi et al.ll2003br 



3 SPECTRAL MODELLING 

The photospheric spectrum obtained on Aug. 17, 2004 (day 
+47, Table 6) is modelled using the gY T VO 14^ parametrized 
spectrum synthesis code llFisheil l2000l : iBaron et al.l l2000l : 
iBranch et alJ l2003l "2004 ) . This code assumes homologous 
expansion of the SN cjccta, i.e. v{r) = r ■ Vph/Rph, where r 
is the radius of a thin shell above the photosphere, v{r) is 
the velocity of the shell, Rph and Vph is the photospheric ra- 
dius and the velocity at the photosphere, respectively. Line 
formation is treated as in the Schuster-Schwarzschild model: 
the photosphere radiates as a blackbody with temperature 
Tbb, and the lines are formed entirely in the envelope above 
the photosphere. The line formation is assumed due to pure 
scattering. This is equivalent to setting the line source func- 
tion as 



Sir) = ^-f 



{^y 



(2) 
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Table 8. Reference line parameters and optical depths of the 
SYNOW model spectrum shown in Fig. 11. The model has been 
computed assuming T;,;, = Ti^xc = 8000 K, Vp^ = 3400 kms~^ 
and o = 6 (see text). The columns give the following data: name 
and ionization state of the element, wavelength (in A), oscillator 
strength, excitation potential (in eV) and the fitted optical depth 
for the reference line. 



Ion 


A (A) 


log(gf) 


x(eV) 


Tref 


HI 


6563 


0.71 


10.21 


200 


Nal 


5890 


0.12 


0.00 


1.1 


Sell 


4247 


0.32 


0.32 


1.0 


Till 


4550 


-0.45 


1.58 


8.0 


Fell 


5018 


-1.40 


2.89 


4.0 


Ball 


4554 


0.19 


0.00 


2.0 



where Iph ~ Bx{Tbb) is the specific intensity at the pho- 
tosphere. The optical depth as a function of radius was ex- 
pressed as a power law: 



iiph 



(3) 



In SYNOW the fit of the spectral lines is controlled 
by adjusting the optical depth of a reference line (a strong 
line in the optical regime) for each element. The depths of 
other lines are computed via the Boltzmann-equation. The 
solution of the radiative tr ansfer equation is performed in 
the Sobolev approximation JKasen et alJl2002f) . 

Table 8 lists the elements identified in the day -1-47 
spectrum, together with their atomic parameters and fitted 
reference line optical depths. The model spectrum and the 
observed one are plotted in Fig. 13. The photospheric tem- 
perature was set as Ttt = 8000 K determined from the slope 
of the red continuum. The velocity at the photosphere was 
chosen as Vph ~ 3400 kms~^, which gives Rph = 922 Mkm 
(1324 Rq) for the radius of the photosphere -1-47 days after 
explosion. The optical depth exponent that gives the best 
description of the observed spectrum was found as a = 6. 

The elements that could be identified unambigously are 
H, Na I, Sc II and Fe II. Ti II is assumed to be responsible 
for the blends below 4800 A, but other elements (e.g. Ca, 
Mg, Ni, Co) may also contribute. Ba II is also present, but 
its contribution is weak and blended by Fe II. Because the 
overall appearance of the spectrum is very similar to that 
of SN 1999em (Fig. 10), the m ore extensive line identifica- 



tion bv lLeonard et al.l ll2002af) is probably also valid in this 
case. Fig. 14 shows the contribution of the elements listed 
in Table 8 to the spectrum of SN 2004dj. 



3.1 The Ha profile 

It is seen in Fig. 13 that Ha cannot be fitted well by 
SYNOW: the absorption is too deep and the emission is 
too low in the model spectrum. This is a usual situation in 
modelling a strong line with a pure scattering source func- 
tion (Branch et al. 20Q3^ . It does not necessarily mean the 
breakdown of the Sobolev approximation, but the source 
function probably differs from that in Eq.l. 

We found that a simple modification of the line source 
function can sufficiently explain the emission content of Ha 




4000 4500 5000 5500 6000 6500 7000 7500 8000 
Wavelength (A) 

Figure 13. The SYNOW model fit (thick line) and the +47 days 
observed spectrum (thin line). The model parameters are listed 
in Table 8. 




4000 4500 5000 5500 6000 6500 7000 7500 8000 
Wavelength (A) 

Figure 14. The contribution of chemical elements to the ob- 
served spectrum. The labels on the right-hand side indicate the 
composition. The adopted parameters are the same as in Table 

8. 



within the Sobolev approximation. Now the source function 
takes the form 



S{r)=K Sn{r) 



(4) 



where So{r) is the function given in Eq.l and _R" is a constant 
in r. In fact, K simply scales the photospheric intensity, 
but the shape of the source function remains the same as 
the pure scattering one. Fig. 15 shows the line profiles com- 
puted with K — 1 and 2. The agreement with the observed 
spectrum is much better when K = 2, i.e. the photospheric 
intensity has an excess with respect to the blackbody in- 
tensity at the wavelength of Ha. Since the photosphere is 
actually a hydrogen recombination front in SNe Up, the ex- 
cess emission can be interpreted as due to the recombination 
of hydrogen to the n — 3 excited state followed by a down- 
ward transition to the n = 2 state, producing spontaneous 
emission at the wavelength of Ha. The exact form of the 
source function in SNe Up atmospheres would deserve fur- 
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Figure 15. The Ha line profile computed with the scaled source 
function {Eq.4). The continuous line corresponds to ii" = 2, while 
the dashed line is for K = 1. The observed line profile is also 
shown for comparison (thick line). 



ther study. iKasen et al.l (|2003) presented a technique for de- 
termining the source function (as weU as the optical depth) 
as a function of r. We leave this interesting problem for fu- 
ture studies. 



4 DISTANCE DETERMINATION 

There are a number of methods for determining SNe lip 
distances. The Expand i ng Photosphere Method (EPM) 
jKirshner fc KwaDlll974l: iHamuv et alJl20oil iLeonard et alJ 



l2002allbt iDessart fc Hillie J 1200531) . 



a variant of the famous 
Baade-Wesselink method, derives the angular diameter from 
both photometry and radial velocities, and their compari- 
son yields the distance. The more sophisticated extension of 
this method is th e Spectral Fitting E xpanding Atmosphere 
Method (SEAM) jBaron et alJll994l) . which uses a custom- 
crafted SN model atmosphere and full treatment of hydro- 
dynamics and radiative transfer to fit model spectra to the 
observed ones. SEAM has been succesfuUy applied to SNe 
Up and other types as well. SEAM stands on more solid 
physical basis than EPM, but requires much more comput- 
ing power. 

T he Standard Candle Method (SCM) jHamuv fc Pintd 
|2002^ relies on the correlation between the absolute magni- 
tude and the radial velocity of the SN ejecta in the middle 
of the plateau phase. Here the distance can be found by 
comparing the observed and the calculated absolute magni- 
tudes, like in the case of Cepheids. The so-called "Plateau- 
Tail Method" (PTM) jNadvozhin|l2003l) compares the SN 
light curve with theoretical models and derive the distance 
via the model parameters. 

In the following we apply EPM and SCM to derive the 
distance to SN 2004dj, and compare it with distances of the 
host galaxy found by other methods. 

4.1 Expanding Photosphere Method 

The methodology of EP M has been recentl y review ed and 
thoroughly discussed by iDessart fc Hillieii ||2005a|). They 



Table 9. Quantities derived during the EPM. The columns con- 
tain the foUowings: JD —2450000 of observation, temperatures (in 
Kelvins), angular sizes (in 10^ km Mpc~^), and values of O/vpi^ 
(in day Mpc"-"^). 



t 


Tbvi 


Tvi 


Obvi 


dvi 


(d/v)BVI 


{e/v)vi 


3223.4 


7047 


10602 


3.347 


1.700 


9.674 


4.912 


3226.6 


6632 


9988 


3.618 


1.967 


10.984 


5.973 


3228.6 


6419 


9795 


3.684 


2.033 


11.568 


6.385 


3229.6 


6412 


9727 


3.687 


2.062 


11.752 


6.571 


3230.4 


6823 


9851 


3.570 


2.043 


11.511 


6.589 


3234.4 


5955 


9653 


3.722 


2.023 


12.841 


6.980 


3236.6 


6079 


9385 


3.729 


2.163 


13.340 


7.740 



concluded that EPM is best used at early phases, when the 
SN ejecta is fully ionized and continuum opacity dominates 
line opacities in the optical. Unfortunately, SN 2004dj could 
not be observed at such early phases, which makes the ap- 
plication of EPM problematic. Nevertheless, we try to use 
EPM for the earliest observed data of SN 2004dj and crit- 
ically compare the results with other distance estimates to 
minimize the systematic errors caused by the possible failure 
of the asssumptions of EPM. 

EPM has been applied in the same way as described 
bv I Vinko et al] (2004) . The angular radius has been derived 
from the UVOIR bolometric light curve (Sect. 2. 1.3) as 



fbol 



eiThT^ff 



(5) 



where fboi is the "observed" bolometric flux, T^ff is an ef- 
fective temperature of the diluted blackbody and (,{T) is 
the dilution factor (determined from SNe Up model atmo- 
spheres) . For each epoch, the effective temperature has been 
estimated by fitting blackbodies to the dereddened BVI and 
VI fiuxes (the R band has been omitted because of the pres- 
ence of Ha emission). 

The dilution factor ^{T) in Eq.5 comes from model at- 
m ospheres. We have adopted the dilution factors derived 
by 
by 

derived new dilution factors that are systematical ly higher 
than those of lEastman. Schmidt fc Kirshnerl il996r . The ef- 
fect of using these dilution factors on the distance will be 
discussed later in this Section. Note that the dilution factors 
depend on the bands used for temperature determination, 
thus, they are slightly different for Tbvi and Tyi- 

Assuming spherically symmetric SN ejecta and homol- 
ogous expansion, for each oberved epoch t, the distance D 
is given by the linear equation 



Eastman. Schmidt fc Kirshneil | l996t) a nd parametrized 
Hamuv et al.l ||200J) . Recentlv IDessart fc Hiuieil l|2005a^ 



t = to + D- 



Vph 



(6) 



where to is the moment of explosion, Vph is the photospheric 
expansion velocity at epoch t. The velocities have been in- 
terpolated from the observed values listed in Table 7. The 
derived parameters relevant for the EPM analysis are listed 
in Table 9. 

The EPM quantities are plotted in Fig. 16. The slope 
and the zero point of the fitted line gives the distance and 
the moment of explosion, respectively (see Eq.6). It is visible 
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Figure 16. Distance determination in EPM via Eq.6. Parameters 
from Tgvi ^re plotted witli filled circles, while those from Tyj 
are represented by open symbols. The fitted lines give the distance 
and the moment of explosion (see text). 



Table 10. Results of EPM analysis. The uncertainties are given 
in parentheses. The quality of the fit (rms) is also shown. 



Bands 



to 
(JD-2450000) 



D 

(Mpc) 



BVI 3187 (3) 

VI 3198 (4) 



3.66 (0.3) 0.77 
4.97 (0.6) 1.14 



that the data from Tbvi and Tyi give quite diflForcnt pa- 
rameters, which is a usual situation in EPM ^Leonard ct al. 
l2002al lbt). 

Which solution is closer to the real one? The concept of 
EPM (diluted blackbody, spherical symmetry, etc.) received 
many criticism recently JLeonard et alJ l2003l : iBaron et alJ 
|2004), leading to a conclusion that EPM is at least "suspi- 
cious" and the results are systematically off from the real 
values. The discussion of all the problems is beyond the 
scope of this paper. Here we would like to show that com- 
bining the EPM results with those of other methods and the 
observations, it may be possible to select an EPM solution 
that is consistent with other results. 

A critical step in EPM is the determination of the effec- 
tive temperature, since it affects the calculated blackbody 
flux as well as the dilution factor. This quantity strongly de- 
pends on the selected wav eleng th bands (Ta b le 9). Recently 
iDessart fc HiUieil (jgOoHi) and iBaron et all (|2TOi) showed 
that such kind of fltted blackbody temperature does not 
represent well the real spectral flux distribution of the SN 
atmosphere. Indeed, comparing the temperatures in Table 9 
and the result of spectral modelling in Sect. 3, it is clear 
that Tbvi is too low, while Tvi is too high with respect to 
the temperature {Ti,ij — 8000 K) that represents the spec- 
trum observed on day -1-47 well. This is illustrated further 
in Fig. 17, where the observed spectrum and the SYNOW 
model computation are plotted together with blackbody flux 
distributions corresponding to Tbvi = 5955 K, Tvi — 9653 
K and the average of these two temperatures. The black- 
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Figure 17. Flux distribution of blackbodies fitted in EPM, 
together with the observed and model spectrum on day -(-47 
(Sect. 3). The drop of flux in the blue region caused by metallic 
line blends (Fig. 14) and the continuum slope in the red can- 
not be fitted by a single blackbody. However, the total flux (but 
not its spectral distriution) may be adequately described by the 
blackbody with Tbvi- 



body curves were scaled according to the angular radius of 
the photosphere found by EPM. 

It is clear that the continuum slope on the red side of the 
observed spectrum is indicative of higher temperature, but a 
blackbody with such a temperature fails to model the sudden 
drop of t he flux below 4000 A. Simi lar fl gures have be e n pre - 
sented bv lDessart fc Hiflieil ll2005al) and lBaron et al.l (|20o3) 
for SN 1999em. There is no doubt that the blackbody as- 
sumption is invalid, even in the very early phases (that were 
actually not observed in the case of SN 2004dj). The cause 
of this disagreement has been already shown in Fig. 14: the 
metallic line blends (especially the Ti II absorption) are re- 
sponsible for the sudden decrease of the flux in the blue. 

Fig. 17 shows that the blackbody with Tbvi has a flux 
distribution that does not give the best flt between 4000 
and 6000 A, but it adeqautely describes both the flux drop 
in the blue and the tail in the red. The two other curves 
that are closer to the real photospheric temperature have 
large excess flux below 4000 A. These curves highly overes- 
timate the photospheric flux in the blue bands, increasing 
the EPM distance. The Tbvi blackbody, however, results 
in a bolometric flux that is closer to the real bolometric lu- 
minosity, even though the spectral distribution is not per- 
fectly matched. Since the purpose of EPM is to estimate 
the distance, rather than to model the spectra, the diluted 
blackbody assumption may be an acceptable compromise, 
if one regards it as the representation of the amount of in- 
tegrated flux of the photosphere, rather than a physically 
correct model of the spectral flux distribution. 

In principle, the application of a proper dilution factor 
should correct the flux for both the temperature mismatch 
between the SN and the blackbody, and the line blanket- 
ing. In this ideal case the same distance could be derived 
from any filter combinations. However, this is true only if 
the model spectrum (that was used to compute the dilution 
factor) perfectly fits the observed one. In reality, the metal- 
lic lines in the blue band can vary significantly from SN to 
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SN, which can be accounted for only with custom-crafted 
model atmospheres. This makes the usage of B band fluxes 
and temperat ure such as Tb v less effective for EPM. 

Recently iDessart fc Hillicr (2005a[l presented a new 
set of dilution factors appropriate for SNe Up, that are 
systematically higher by 10 - 20 % than those given by 
lEastman. Schmidt fc Kirshnerl (Il99dl . In a follow-up pa- 
per IPessart fc Hillien |2005b) showed that the application 
of their new dilution factors for SN 1999em increases its 
EPM-distance, bringing it into good agreement with the 
Cepheid-distance of the host galaxy, thus, eliminating most 
of the systematic errors of EPM found by Leonard c^lJ 
('2003). If one applies the dilution factors of lDessart fc Hillien 
(L2005a1 for the present data of SN 2004dj, its EPM-distance 
will increase as Dbvi ~ 4.3 ± 0.3 and Dyi = 6.2 ± 0.7 
Mpc. As in the case of SN 1999em, these distances are also 
systematically higher than those that were derived from 
the Eastman et al. dilution factors. The problem is that 
for SN 2004dj the EPM-distances are all longer than the 
values found by other methods (see bel ow). In this case 
the application of the dilution factors oflD essart fc Hillieil 
( 2005a) even increases the discrepancy between EPM and 
the other methods. It might be possible that t he SNe Up 
models bv Eastman. Schmidt fc Kirshnen il996l) represents 
S N 2004di better than SN 19 99em, while the models apphed 
bv'Pess art fc HiUieil ll2005al) are better for SN 1999em. On 
the other hand, regarding the similarity of the optical spec- 
tra of these two SNe (Fig. 10), it is hard to believe that the 
atmospheres are so different. Detailed modelling of the spec- 
tra of SN 2004dj may help to resolve this discrepancy. 

Concluding the EPM analysis, we adopt the solution 
based on the Eastman et al. dilution factors and Tbvi (Ta- 
ble 10), namely to = 2453187±3 and D = 3.66±0.3 Mpc as 
the moment of explosion and distance, respectively. In the 
foUowings we show that this distance is in agreement with 
other distance estimates for the host galaxy and SN 2004dj 
itself. 



4.2 Standard Candle Method 

The fo rmulae for the SCM have been selected from iHamuvl 
J2005l) as follows: 



Table 11. Comparison of distance estimates for SN 2004dj and 
NGC 2403. 



^^50-^^+ 6.564 log( 



5000' 



51og(i/o--D) - 1.478 (7) 



t'so 



/50-A/ + 5.8691og(— ^) = 51og(iyo-D)- 1.926 (8) 

OUUU 

where the observed quantities (Vso, ^50, «5o) are the V and 
I magnitudes and the radial velocity (in km/s) 50 days past 
explosion (i.e. in the middle of the plateau phase). Ideally, 
the distances obtained from the V and I data should be 
the same, therefore this method can also be used to test the 
consistency of the applied reddening correction. 

From Table 3 and Table 6 we have chosen V50 = 12.04± 
0.03, /so = 11.40 ±0.03 and W50 = 3250 ±250 km/s. The ab- 
sorption corrections were computed from E{B ~ V) — 0.07 
mag and the galactic reddening law, as before. For the 
Hubble constant Ho = 73 km/s/Mpc has been adopted 
JRiess et al.ll2005l) . This value is now based on both SNe la 
and Cepheid variables, and brings the previously conflict- 
ing distance scales into agreement. From these data we ob- 
tained Dv = 3.47 ± 0.2 Mpc and Di = 3.52 ± 0.2 Mpc 



Method 



D (Mpc) 



Reference 



TuUy-Fisher 
Cepheid P-L 
EPM 
SCM 



3.50 ±0.31 
3.22 ±0.15 
3.66 ±0.30 
3.50 ±0.20 



LEDA: Russell (2002) 
Freedman et al. ('20011 



present paper 
present paper 



average 



3.47 ±0.29 



from the V and / calibration, respectively. The uncertain- 
ties are purely from the observational errors. These distances 
agree nicely with each other and with the EPM distance 
{Depai ~ 3.66 ± 0.3 Mpc) within the errors, strengthening 
the credibility of both of the distance and reddening esti- 
mates. 



4.3 Comparison with other distance 
determinations 

The distances computed above are all based on the photo- 
metric and spectroscopic behaviour of the SN ejecta during 
the plateau phase. The overall unce rtainty of such distance s 
is about 20 % in the case of SCM llHamuv fc Pintdl2002l) 
while it can be as large as 60 % for EPM JLeonard et alJ 
I2OO3). Thus, the SN-distances should be compared with 
other distance estimates based on other methods, in order 
to test the possible systematic errors. 

Since the host galaxy, NGC 2403, a member of the M81 
group, is a nearby, well-studied object, its distance has been 
estimated from many other methods. These are listed in 
Table 11 together with the SN-based distances. 

The Cepheid distance JFreedman et al.ll200ll) was usu- 
ally applied in previous studies of SN 2004dj, but it is worth 
noting that this value is based on o nly I photom etry with- 
out direct reddening determination. iFreedman et al. (200J) 
assumed Aj — 0.2 mag as the 7-band absorption during the 
distance measurement. From Table 11 it is visible that the 
SN-distances are in agreement with the recent TuUy-Fisher 
distance, and differ less than 2a from the Cepheid distance. 
This is not a significant difference, if one takes into account 
the additional uncertainty of the Cepheid distance due to 
the lack of reddening information. We conclude that the av- 
erage distance of NGC 2403, 

D = 3.47 ±0.29 Mpc 

is consistent with all available distance estimates within the 
errors. 



5 PHYSICAL PARAMETERS 

5.1 Nickel mass 

During the tail phase the light curve is thought to be pow- 
ered by radioactive decay of ^''Co - ^''Fe. Since ^''Co is the 
daughter nucleus of ^^Ni, the luminosity is proportional to 
the initial nickel mass produced by the explosion. 

In order to estimate the nickel mass of SN 2004dj we 
compared the absolute bolometric light curve (calculated 
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Figure 18. Comparison of absolute bolomotric magnitudes of SN 
2004dj (filled circles) with model computations (lines) during the 
early tail phase. Labels indicate the Ni-masses (in Mq). 



Table 12. The adopted magnitudes for Sandage 96 (see Section 
2.1 and Table 2). 

Filter A (A) Mag. Reference 



u 


3600 


17.68 






B 


4400 


18.26 


present paper 




V 


5500 


17.85 


present paper 




R 


7000 


17.53 


present paper 




I 


9000 


17.06 


present paper 




J 


12500 


16.19 




ri997'l 


H 


16500 


15.54 


Sckrutskie et al. 


(1997) 


K 


22000 


15.42 


Sckrutski^^^L 


(19971 



is MNi = 0.02 ± 0.01 Mq. This is sli fihtly below the 
nickel mass of SN 1999em and SN 1999gi (lElmhamdi et alJ 
|2QQ2a|), which is an expected result based on the light curves 
(Sect.2.1.1). 



from the "observed" UVOIR bolometric light curve dis- 
cussed in Sect. 2. 1.3, and the distance estimated in Sect. 4) 
wi th the predict i on of the simple light curve model described 
bv lVinko et alJ l|200^. This model calculates the luminos- 
ity in a homologously expanding SN ejecta, assuming ra- 
dioactive decay, full gamma-ray and positron trapping and 
prompt thermalization of gamma photons and positrons. 
The density structure of the ejecta has been given as a 
power- law with exponent —7. This value is derived from the 
radius dependence of the optical depth r(r) = To{r/Rph)~^ 
(found in Sect. 4) by computing p{r) — — k,^^ (dr / dr) and 
assuming spatially constant k. The ejected mass has been set 
as Mej = 2OM0 (see Sect. 5.2). The calculated light curves 
were mostly sensitive to the initial nickel mass, which was 
treated as a free parameter. Fig. 18 shows the light curves 
(thin lines) calculated with Mni = 0.01, 0.02 and 0.03 Mq. 
The observational data are also plotted for comparison as 
filled circles. Although during the observations SN 2004dj 
did not enter into the radioactive tail phase completely (the 
slope of the observed light curve differs from the calculated 
one given by the Co-decay), it is seen that the luminos- 



0.02 ± 0.01 M( 



SL, 



ity level can be explained with Mjv 

This nickel mass is the same as derived by IChueai et al.l 
l|2005|) from the comparison of the absolute V light curves 
of SN 2004dj and SN 1987A. Based on the steepness of the 
V light curve measured after 100 davs iKotak et aU ll2005l) 
got MNi ~ 0.022 Mq. 

Alternatively, the nickel mass can be determined us- 
ing semi-empirical corre lations, such as those presented by 
lElmhamdi et al.1 1|2003^. Using their Eq.2 one can compute 
the nickel mass from the absolute V magnitude 35 days 
before the light curve inflection point. The moment of in- 
flection point was found as Ti = 2453297 JD (Sect.2.1.1), 
and the absolute magnitude 35 days before that is Mv {Ti — 
35) = —15.52 ± 0.1 mag. The nickel mass is then MMi = 
0.022 ± 0.002 Mq, in perfect agreement with the value de- 
rived above. From the steepness of t h e ligh t curve dur- 
ing the transition phase IChueai et al.l (I2OO3) got Mjvi = 
0.013 ± 0.004Mo. It is concluded that aU the methods ap- 
plied here agree that the estimated nickel mass of SN 2004dj 



5.2 Progenitor mass 

Up to now the progenitors of Type II SNe have been identi- 
fied directly in a number of cases (see Sect.l for references). 
Surprisingly, these turned out to be relatively less massive 
stars, some of them being close to the theoretical limit of 
core collapse (~ 8Mq). However, in most cases the progen- 
itor could be detected only in a single bandpass, preventing 
the identification of its colour and evolutionary state. 

SN 2004dj is particularly interesting in this respect, be- 
cause its progenitor was a member of a compact cluster 
(Sandage 96) in NGC 2403. The age of a cluster can be 
determined more precisely than that of a single star, thus, 
the evolutionary state of the progenitor could be secured. 
Two groups have reported their observations of S96, made 
some years before the explosion JMaiz-Apellaniz et alJl2004l : 
IWang et alJl2005h and their analysis of the cluster SED (see 
Sect.l). 

We have re-analysed the SED of S96 based on our own 
BVRI ph otometry (Table 2) supplemented by data from the 
literature iLarsenlll999l : ISckrutskie et al.lll997l) . Our optical 
photometry is in ver y good a greeme n t with the observations 
published by iMaiz-A pcUani z et al.l ll2004r) . except for the 
B — V colour ind e x, wh ere our value is closer to the one given 
bv lChueai et alJ ll2005ll . The a dopted magnitudes of S96 are 
shown in Table 12. Note that iMaiz-Apellaniz et al.1 (|20o3) 
reported a Gunn- Johnson color index u — B — 0.41 ± 0.06 
that d iffers from the U ~ B — —0.51 value given bv iLarsenl 
il999h in the Johnson system. We have preferred the lat- 
ter value, since it is in the same photometric system as the 
other optical fluxes, but the U magnitude seems to be more 
uncertain than the fluxes in other bands. It may affect the 
derived cluster age, since it is connected with the amplitude 
of the Balmer-jump that increases with age. 

We have selected two recent grids for calculating 
the model SEDs: the w idely used model sequence by 
iBruzual fc CharlotI ||200J ) (BC03), based on Padov a evol u- 



tionar y tracks (that was also applied for S96 bviWang et al. 
120051)), and a more recent one published bv I Jimenez et al, 
1120041) (J04), containing new stellar interior models, evolu- 
tionary tracks and improved treatment of mass loss. For the 
cluster initial mass function, a Salpeter IMF was adopted. It 
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Figure 19. Results of SED fitting to the observed data of S96. The applied grids and the fitting parameters are indicated on each panel. 
See Table 13 for more information. 



is known that the Salpeter function does not describe well 
the shape of the IMF, especially in the low mass regime. 
The selection of the Salpeter IMF was motivated mainly by 
the fact that the two previous studies for S96 also used this 
function, and we wanted to make the comparison of different 
solutions easier. For the same reason, a single starburst was 
assumed as the star formation history of S96. 

The fitting has been computed via a simple x^ mini- 
mization. The observed magnitudes (Table 12) were trans- 
formed into monoc hromatic fluxe s (in e rg/s/cm^/A) using 
the calibration by iHamuv et alJ (|200l|) (see Sect. 2.1.3), 
dereddened, and compared with the predictions of the spec- 
tral synthesis model grids. The fitted parameters were the 
cluster age Td, the cluster mass Md and the overall cluster 
reddening E{B — V)ci- The metallicity of the cluster was 
kept fixed during the minimization. It should be noted that 
E[B — V)ci niay be quite different from the reddening de- 
rived for SN 2004dj, because the compact cluster may con- 
tain a significant amount of intracluster gas and dust that 
is situated behind the SN. On the other hand, these pa- 
rameters are correlated, and their simultaneous fitting may 
lead to spurious solutions. This is especially true for the red- 
dening and the age, because a higher reddening correction 
would result in a bluer SED, indicating a lower age. 

The resulting parameters corresponding to the min- 
imum of x^ ^re listed in Table 13 and plotted in 
Fig. 19. Because the metallicity of NGC 2403 at the 
position of S96 is probably sub s olar ([0/H1 = — 0.24, 
IPilvuein. Vilchez fc Continil'^l2004^ : IWang et ai](l2005^ ). we 



have used the models with Z = 0.02 and Z — 0.004 metal- 
licities from both grids. Comparing these results with the 
published ones (shown in the last three rows of Table 13), 
it is visible that our analysis resulted in the same kinds of 
solutions as the previous studies: a "young" population with 
Trf = 8 Myr and an "old" population with 20 < T^i < 30 
Myr. The "old" solutions are very similar to the ones given 
in the previous studies. The "young" solution, however, has 
considerably less age. It is interesting that this solution 
comes out using both model grids, regardless of metallicity 
(but with different reddenings and cluster masses, of course). 
Based on the x^ values, the "young" solution from the BC03 
grid gives the best description of the observed data. How- 
ever, we agree with lWang et alJ i2005fl that if the observed 
SED is undersampled, such as the one presented here, it 
may be misleading to accept or reject solutions based on x^ 
values only. Thus, we conclude that, at present, there are 
three probable solutions for the SED fitting of S96: the old 
one corresponds to 20 < T^i < 30 Myr, the young one with 
Tci ~ 8 Myr, and a solution with age between these two, 
Tci = 14 Myr. Each solutions have their own combinations 
of cluster reddening and total mass, between the values of 
0.04 < E{B -V)= 0.35 and 30 < M^ < 140 ■ 10^ Mq. 
Due to the correlations between these parameters, it is not 
possible to assign a unique solution for all these quantities 
simultaneously. 

Regarding the potential progenitor of SN 2004dj, it 
is seen in Table 13 th at the SN ma s s (de termined from 
the Padova isochrones, iGirardi et alJ il996f) l ranges from 
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Table 13. Results of SED fitting for S96. The columns contain the foUowings: Model sequence, metallicity, cluster age (in Myr), cluster 
mass (in 10'^ Mq), cluster reddening (in mag), predicted SN mass for the giv en cluster age (in M(r) , estim ated from Padova isoch rones), 
reduced x^ and reference. The last three rows list the solutions published bv lMaiz-Apellaniz et al.1 i2004h and lWang et al.1 J2n' 
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~ 10 Mq to ~ 25 Mq, depending on age. The uncer- 
tainty of the age determin ed from SED fitting is at least 
10 % l|,Timenez et al J 1200 J . implying an error of the SN- 
mass of about ±4 Mq for the young- and ±1 Mq for the 
old solution. It is interesting that the young (8 Myr) solu- 
tion results in significantly higher mass than all the SNe 
Up progenitor masses detected so far. The highest mass 
of a detected SN Type II progenitor is that of SN 1999ev 
JMaund fc S martt 2005) with 15 Mq < Mzams < 20 M©. 
This is similar to the pro genitor mass of SN 2004et, be- 
ing Mzams = 15^2 Mq ilLi et alJl2005a^ and SN 1999gi 
with Mzams = 15t| Mn. iJLeonard et alJ l2002bh. Con- 
trary to these, the progenitors of SN 2003gd (Van Dvk et al.l 
[2fl03; Smartt et al. 2004) and SN 2005cs ( Li et al. 20051a; 
iMaund. Smartt fc Danzigerl l200,4) . also detected on pre- 
explosion images, are found within the mass range of 
Mzams = 7—9 Mq. These two recent discoveries have led to 
a conclusion that perhaps all Type Up SNe emerge from rela- 
tively less massive progenitors, so that Mzams ~ 8 — 15 Mq. 
This is supported by the non-detection of the progenitors o f 
SN 1999em, SN 1999gi and SN 2001du iSmartt et alJl2003l) 
giving Mzams < 15 Mq . 

However, the possible masses of the progenitor of 
SN 2004dj determined above do not fully support this sce- 
nario. If the Tci = 8 Myr solution of S96 is indeed valid, then 
the progenitor mass of SN 2004dj must have been higher 
than 20 A'Iq. This is significantly higher than the masses 
found for the other progenitors, even for SN 1999ev (see 
above) . Since the other two solutions for S96 give progenitor 
masses that are consistent with the mass range of the oth- 
ers, it cannot be stated that SN 2004dj had definitely a more 
massive progenitor, but it cannot be ruled out at present. In 
principle, the cluster masses in Table 13 imply that there 
might be at least a few stars with M > 20 Mq when 
SN 2004dj exploded. From the integration of the Salpeter 
IMF, the number of stars with M > 20 Mq is N{M > 
20) « 2 • (M/W^Mq), which gives N{M > 20) = 58 for 
Mci = 29 • 10^ Mq. The actual number of massive stars 
in S96 may, of course, differ from this value, but statisti- 
cally it is possible that at least a few such stars are indeed 
there. [Maiz-Apellaniz et al. ( 2004i) estimated that S96 con- 
tains about 12 red and perhaps 2-3 blue supergiants, based 
on their SED parameters. The spectrophotometric observa- 
tions of the remainder of S96 after the SN faded away would 
help to clarify the picture. 



Comparing the observations of the CO vibrational 
bands in the mid-IR with theoretical explo sion models 
combined with radiative transfer calculations, iKotak et alJ 
(2005) presented evidence for a red supergiant with ~ 15 Mq 
as the likely progenitor of SN 2004dj. Again, this is at the 
upper end of the detected progenitor masses for Type II SNe, 
close to the prediction of th e ~ 14 Myr cluster solution by 
' Maiz- Apellaniz et al.' 1*2004'). 



|Zhang et al. (2005) estimated the parame ters of the 
are-su pernova star using the formulae given bv iNadvozhinI 
i2003r . These formulae relate the length of the plateau 
Atp, the middle-plateau absolute magnitude Mv and ex- 
pansion velocity at the photosphere Vph to the explosion 
ener gy (Eexv), ejected m ass {Mej) and progenitor radius 
(ffl. lzhanget alJ ^M^ derived Atp = 80 ± 21 days, 
Mv = -16.55 ± 0.35 mag and Vph = 3933 ± 189 kms"\ 
and they got E^^p = 0.751!^^^ x 10^^ erg, Mej = Wtl Mq 
and R = 282^122 ^0- From the more detailed analysis 
presented in this paper, we estimate these parameters as 
Atp = 100 ± 20 days (Sect. 2.1.1, allowi ng 6t » 10 days f or 
the peak duration before the plateau, see lNadvozhinI i2003i) 1 . 
Mv = -15.88 ± 0.18 mag and Vph = 3250 ± 250 kms"^ 
(Sect. 4.2 and 4.3). These parameters result in Eexp ~ 



0.86: 



i X 10''^ erg, Mej « lOI^o Mq and R « 155: 



R. 



It is interesting that the ejected mass is, again, close to 
~ 20 Mq, although with high uncertainty. It is concluded 
that the presently available information suggest that the 
progenitor mass of SN 2004dj is Mprog > 15 M©. The SED 
fitting of S96 may suggest a cluster age as low as 8 Myr, in- 
dicating Mprog > 2QMq. This massive progenitor may also 
be suspected from theoretical SN models. 



6 CONCLUSIONS 

The conclusions of this paper are summarized as follows: 

• New BVRI photometric observations of SN 2004dj are 
presented. The date of explosion is estimated as June 30, 
2004 (JD 2453187±20), about 1 month before discovery. The 
plateau phase lasted about ~ 110 ±20 days after explosion. 
The most probable value of the re ddening toward SN 2004dj 
is E{B -V) = 0.07 ± 0.1 mag JCuenther fc Klosdl2004^ . 
The progenitor cluster S96 is also detected on pre-explosion 
frames. 
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• The new, plateau-phase optical spectra reveal great 
similarity between SN 2004dj and SN 1999em, a typical SN 
lip. Radial velocities from metallic lines and Ha are pre- 
sented. The spectral features are successfully modeled with 
the spectrum synthesis code SYNOW with Tbb = 8000 K as 
the photospheric temperature. The observed P Cyg profile 
of Ha can be better described by an increased emission of 
the photosphere relative to the blackbody intensity at the 
rest wavelength of Ha. A single nebular spectrum c onfirms 
the Ha asymmetry reported bv lChueai et al.l (l2005h . 

• The distance to SN 2004dj is estimated from two meth- 
ods applicable for SNe Up. The EPM distance is Depm ~ 
3.66±0.3 Mpc, while the SCM gives Dscm = 3.50±0.2 Mpc. 
These are in good agreement with other distance estimates 
to the host galaxy NGC 2403, based on Cepheids and TF- 
relation, although the Cepheid-distance is slightly less than 
the SN-distances. The averaged distance is D = 3.47 ± 0.3 
Mpc, which is consistent with all available information. 

• Using the updated distance, the nickel mass of 
SN 2004dj is calculated as Afiv, = 0.02 ± 0.01 M© from 
the shape of the tail light curve. The progenitor mass has 
been estimated from fitting theoretical SEDs to the ob- 
servations of S96. The so lutio ns are similar to those of 
iMaiz-Apellaniz et alJ (|200i) and lWang et alJ (120051') . except 
that a "young" solution corresponding to Td = 8 Myr is 
also possible. Such a young cluster age would mean that the 
mass of the progenitor was Af > 20 Mq. The SN-masses 
of the other solutions are less than 15 Mq , being in better 
agreement with other SNe Up progenitors. 
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